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I discuss CLEO-c opportunities in Dalitz plot analyzes with the data samples available 
now and projected by the end of CESR-c run. Using 281 pb~^ of e+e" collisions at mass 
of i/i(3770) we present results of the Dalitz plot analysis of D+ — > 7r~7r"'"7r+. Using the 
CLEO-c and CLEO III samples of 5 pb~^, accrued at mass of ip{2S), we study three 
body decays of XcJ produced in the radiative decay ip{2S) — > 7XcJi where J=0,l,2. 
A clear signal from at least one of Xaj is found in eight final states: 7r+7r~r7, TT'^n~r]' , 
K+K-n°, K^K-K, r]K+K- ' , TT^pp, r]pp, and AK^p. For these modes we measured or set 
an upper limit on the branching fraction. A Dalitz plot analysis is performed on three 
modes Xci — > n^n~ri, K~ , and K^Ktt. 
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1. The CLEO-c opportunities in Dalitz plot analyzes 

Currently the CLEO collaboration possesses multiple samples of data which can be 
used for analysis of D, Dg, Xcj, J/ip meson decays. 

• A sample of e+e" — > tp(377Q) — > DD events with integrated luminosity of 
281 pb-\ which is equivalent to - 10^ of D+D-, and 0.8 • lO*' of D^'W pairs. 
We find about 20 three-body decay modes of charged and neutral D mesons to 
the combinations of tt, K, and ?/ in the final state that can be used for the Dalitz 
plot analysis. Running at DD production threshold gives many advantages. There 
is no extra-energy for production of additional fragmentation particles which keeps 
the multiplicity low and gives very clean events. In some cases, where combinatoric 
background is small, we use inclusive reconstruction of a single D meson. In exclu- 
sive modes one D meson of DD pair can be used as a flavor or CP tag^, and the 

*Presented at Charm-2006 conference, June 5-7, IHEP, Beijing, China. 
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recoiling D studied. 

• A sample of e+e" -> 0*17^ events at \/s = 4170 MeV with C ~200 pb^^, or 
2 • 10^ of D*Ds (charge conjugated modes are implied throughout this paper). This 
sample has been collected during Spring of this year and its processing is still in 
progress. We are going to use mesons for the Dalitz plot analysis. In particular, 
for the decay mode I?+ (/)7r+ , often used as a reference, an interference with other 
modes can be studied in the amplitude analysis of the Df K^K~tt^ decay. 

• A sample of e+e" "0(25') events with C ~5 pb~^, or 3 million produced ?/'(2S')'s 
which can be used to study J/ip or Xcj decays in the processes 

l> ip{2S) nnJ/ip with B ~ 50% for tt+tt" and T:°Tr° together, 
O tpi'^SI) ^XcJ with B ~ 9% for each spin J = 0, 1, 2 state. 

The TTTT pair in the first case and the 7 in the second can be used as an additional 
signal tag for these processes. A partial wave analysis can be applied in order to 
study decay of particle with non-zero spin, but as a first step, with this small 
sample in hand, we perform Dalitz plot analyzes with a special treatment of angular 
distributions. 

The current CLEO-c run plan has running on the ^{2S) completed during the 
fall of 2006 accumulating ten times more statistics of ip{2S). By the end of the 
CLEO-c program (April 2008) we are going to gather samples of DD three and 
D*Ds four times larger than we have at present. This data size and cleanliness is 
encouraging us to proceed with decay substructure study in 3-body decay modes 
using the Dalitz plot technique El. At this conference we present preliminary results 
of a few Dalitz plot analyzes of the decays £>+ 7r~7r+7r+ and xa 7r+7r~ry, 
K+K^TT^, and K'^Kn. 

2. Dalitz plot analysis of £)+ — > 7r~7r+7r+ 

A Dalitz plot analysis of 7r~7r'^7r+ has previously been done by E791 ^ 

and FOCUS The preliminary analysis described here is from CLEO-c and 
represents the first time we have done the same Dalitz plot analysis as the fixed 
target experiments. The decay is selected with cuts on the beam constrained mass 
of three charged tracks consistent with pions, rriBc = y/^'beam — P^('?r^7r+7r+), and 
the difference of their energy from the beam energy, AE — E{tt^tt^tt^) — Eb^am- 
A sample of 6991 events is selected with a signal to noise of about two to one. The 
E791 and FOCUS samples are of similar size and cleanliness. 

The Dalitz plot is symmetric under the interchange of like-sign pions thus we 
do the fit in the two dimensions of high versus low unlike-sign pion mass, as shown 
in Fig. njleft). There is a large contribution from £>+ ^^s''^^ which because of 
the long Ka, lifetime should not interfere with other submodes. We do not consider 
events with m(7r+7r~) within ten standard deviations in mass resolution of the 
known K^^ mass on the Dalitz plot. This leaves 4068 events which are shown in the 
projection plots of Fig. ^ An unbinned maximum likelihood fit is used along with 
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Fig. 1. Dalitz plot and projections for — > tt tt+tt 



other methods as described in Ref.l^. The efficiency across the Dahtz plot is modeled 
with simulated events that are fit to a two-dimensional second order polynomial 
and threshold functions to account the fall off at the corners. Backgrounds are 
taken from msc a-nd AE sidebands and extra resonance contributions are allowed 
from mis-measured K^, p, and /2(1270) decays. A total of 13 different resonances 
are considered. Parameters describing these resonances are taken from previous 
experiments. Only contributions with an amplitude significant at more than three 
standard deviations are said to be observed, and others are limited. Contributions 
that are not significant are not included in the decay model used for the result. 

Contributions from p^Tr, /o(980)7r, and /2(1270)7r are clearly visible. Tabled 
shows the preliminary amplitudes, phases, and fit fractions measured by CLEO 
compared with the results of the E791 analysis mentioned above. There is broad 

Table 1. Summary of results for D+ — + 7r~7r+7r+. 



Mode 


Amplitude 


Phase (°) 


Fit Fraction {%) 


E791: Fit Fraction (%) 


p(770)7r+ 


1 (fixed) 


0{fixed) 


20.0±2.3±0.9 


33.6±3.9 


/o(980)7r+ 


1.4±0.2±0.2 


12±10±5 


4.1±0.9±0.3 


6.2±1.4 


/2(1270)7r+ 


2.1±0.2±0.1 


-123±6±3 


18.2±2.6±0.7 


19.4±2.5 


/o(1370)7r+ 


1.3±0.4±0.2 


-21±15±14 


2.6±1.8±0.6 


2.3±1.7 


/o(1500)7r+ 


1.1±0.3±0.2 


-44±13±16 


3.4±1.0±0.8 




(T7r+ 


3.7±0.3±0.2 


-3±4±2 


41.8±1.4±2.5 


46.3±9.2 


Non-resonant 






<3.5 


7.8±6.6 


p(1450)7r+ 






<2.4 


0.7±0.8 






90.2 


116 



agreement between the two results, including the observation of a air contribution. 
In an alternative fit using the same decay model as E791 the agreement is only 
slightly better, but the fit is much less likely to model our data than the model shown 
in the table which does not include non-resonant and p(1450)7r+ contributions, but 
does include a /o(1500)7r+ contribution. Models without a an contribution do not 
agree well with the data. 

This CLEO analysis is preliminary, and we plan to consider a generalized model 



7, 2008 8:29 WSPC/INSTRUCTION FILE text 



4 

of TT+TT^ S-wave interactions to model a and /o contributions such as the K-matrix 
which is used in the FOCUS analysis mentioned above. 

3. Three-body decays of XcJ 

Decays of the Xcj states are not well studied both experimentally and theoretically. 
Assuming that XcJ are the ^Pj cc-bound states one would expect that Xco and 
with J^'-^ quantum numbers 0"*"+ and 2+"*" decay to the light quarks via two gluons^. 
Measurement of any possible XcJ hadronic decays provides valuable information on 
possible glueball dynamics. 

We study the decay modes of Xcj using the process e+e~ '0(25') IXcJ with 
3 million ip{2S) produced. We search for signals in the XcJ h'^h^h^ three-body 
decays, where /i° stands for 7r°, K'^, rj, rj' or A, and /i+ is 7r+, or p. 

Our basic technique is an exclusive analysis. A photon candidate is combined 
with three hadrons and the 4-momentum sum constrained to the known beam 
energy and small momentum caused by the beam crossing angle. We cut on the x^ 
of this fit, which has four degrees of freedom, as it strongly discriminates between 
background and signal. Efficiencies and backgrounds are studied in a GEANT-based 
simulation of the detector's full response. Our simulated sample is roughly ten times 
our data sample. The simulation is generated with a 1 + A cos^ 9 distribution in 
cos 9, where 9 is the radiated photon angle relative to the positron beam axis. A El 
transition, as expected for ^{2S) — > 7Xcj, imphes A = 1, —1/3, +1/13 for J = 0, 1, 2 
particles. The differences of efficiencies due to various 9 distributions are negligible 
as we accept transition photons down to our detection limit. 

We use standard CLEO algorithms for tracks reconstruction in the tracking 
system and photons in Csl calorimeter. Combined dE/dx and RICH, E/p are 
exploited for particle identification and background suppression, tt'^'s are formed 
from two photon candidates of the tt'' 77 decay. We use rj —> 77 and rj tt+tt^tt^ 
decays to reconstruct 77 mesons. The 77' candidates are formed from 77' — s- 7]tt^'k~ 
and and 77' 7p (with E^ > 200 MeV, and tt+tt" mass to be within 100 MeV/c^ 
of the p mass) combinations. In most cases mass constrained fits are applied with 
appropriate cut on value. Kg tt+tt" and A pTT~ candidates are formed 
from good quality tracks that are constrained to come from a common vertex. The 
Kg flight path is required to be greater than 5 mm and the A flight path greater 
than 3 mm. The mass cut around the Kg mass is ±10 MeV/c^, and around the A 
mass ±5 MeV/c^. Events with only the correct number of selected tracks for the 
decay mode being considered are accepted. 

The efficiencies averaged over the CLEO III and CLEO-c data sets for each 
mode, including the branching fractions of 77 — > 77, 77 tt+tt^tt" and 77' t^tt+tt", 
are in the range 8% to 30%. 

Figure 121 shows the mass distributions with fits for the eight Xcj decay modes. 
Signals are evident in all three Xcj states, but not in all the modes. Backgrounds are 
small. The mass distributions are fit to three signal shapes, Breit-Wigners convolved 
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Fig. 2. Mass distribution for candidate Xcj ~* tt+tt rj, K'^K rj, pprj, tt+tt rj' , pjm*^ , K^pK, 
K+K-ir°, and K^K-tt+. 



with Gaussian detector resolutions (with <t from 4.3 to 7 MeV/c^, depending on 
mode), and a hnear background. The Xcj masses and intrinsic widths are fixed at 
the values from Ref. 

We consider various sources of systematic uncertainties on the yields: allowing 
the Xcj masses and intrinsic widths to float in the fit; change the shape of back- 
ground; break up the sample into CLEO III and CLEO-c data sets; uncertainty 
on the efficiency of photon and track reconstruction; uncertainty due to the cut on 
the in mass and event vertex constrained fits. Based on the results of the Dalitz 
plot analyzes we correct the efficiency in the Xci ^ tt+tt^t], Xci ^ K^K^ir'^, and 
Xci — > KgK^TT^ modes by a relative up to 6% to account for the change in the 
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efficiency caused by the deviation from a uniform phase space distribution of de- 
cay products to what we actually observe. To calculate Xcj branching fractions, we 
use previous CLEO measurements 1^ for the '0(25') ^Xcj branching fractions of 
9.33± .14± .61%, 9.07± .11 ± .54% and 9.22± .11 ± .46% for J = 0, 1, 2 respectively. 
Preliminary results for the three body branching fractions are shown in Table |21 
Where the yields do not show clear signals we calculate 90% confidence level upper 
limits accounting for statistical and systematic uncertainties. 



Table 2. Preliminary branching fractions in %. Uncertainties are statistical, systematic, and a sepa- 
rate systematic due to uncertainties in the •il){2S) branching fractions. Limits are at the 90% confidence 
level. 



Mode 


XcO 


Xcl 


Xc2 






< 0.021 




0.52 ± .03 ± .03 ± .03 


0.051 ± .011 ± .004 ± .003 


K+K-q 




< 0.024 




0.034 ± .010 ± .003 ± .002 


< 0.033 


ppv 


0.038 ± 


.010±.003± 


.02 


< 0.015 


.019 ± .007 ± .002 ± .002 






< 0.038 




0.24 ±.03 ±.02 ±.02 


< 0.053 






< 0.006 




0.200 ±.015 ±.018 ±.014 


0.032 ± .007 ± .002 ± .002 


K°K--K+ 




< 0.010 




0.84 ± .05 ± .06 ± .05 


0.15 ±.02 ±.01 ±.01 




0.059 ± .010 ± .006 ± 


004 


0.014 ± .005 ± .001 ± .001 


0.045 ± .007 ± .004 ± .003 


K+pA 


0.114 ± 


016 ± .009 ± 


007 


0.034 ± .009 ± .003 ± .002 


0.088 ±.014 ±.07 ±.006 



4. Dalitz plot analysis of Xci — *■ tt+tt t], K~^K tt", and K^Ktv 

We choose the three high statistics signals Xci — > Tr+vr^r/, Xci —* K^K^tt^ , and 
Xcl KgK^TT^ with 228, 137, and 234 events respectively for Dalitz plot analyzes 
to study resonance substructure. Events were selected within 10 MeV, roughly two 
standard deviations, of the signal peak. When fitting the data Dalitz plot the small 
contributions from backgrounds are neglected. We use an isobar model to describe 
resonance contributions taking into account spin and width dependent effects. Nar- 
row resonances are described with a Breit-Wigner amplitude with the resonance 
parameters taken from previous experiments ■ . We use the Fl atte line-shape for 
ao(980) with parameters from the Crystal Barrel Collaboration For low mass 
7r^7r~((T) and K'k{k) S-wave c ontr ibutions we choose a simple description which is 
adequate for our small sample 

We are examining the e+e" tpi'^S) — > 7Xci process. In such a decay the Xci 
should be polarized. For our small sample we use angular distributions from^l^l cal- 
culated for non-polarized Xci decay. We have tested different angular distributions 
and include the variations as a systematic uncertainty. 

Fig. 13 shows the Dalitz plot and three projections for Xci TT^n^ri. There 
are clear contributions from ao(980)^7r^ and /2(1270)77 intermediate states, and 
significant accumulation of events at low tt+tt^ mass. Isospin symmetry predicts 
that amplitudes and strong phases of both charge conjugated states should be equal. 
The overall amplitude normalization and one phase are arbitrary parameters and 
we set aa(,(980)+ = aao(980)- = 1, 0ao(98O)+ = 4>ao{980)- 0. AU Other fit components 
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Fig. 3. Dalitz plot and projections on the two mass squared combinations for Xci tt+tt rj. 



are defined with respect to this choice for ao(980). 

Our initial fit to this mode includes only ao(980)*7rT and /2(1270)r7 contribu- 
tions, but has a low probability of describing the data, 0.13%, due to the accumu- 
lation of events at low tt+tt" mass. We find that only the a-r] describes well the 
low TT+TT" mass spectrum, and finally describe the Dalitz plot with ao(980)±7rT, 



/2 (1270)77, and arj contributions. Table 13 gives the preliminary results of this fit 
which has a probability to match the data of 58%. Systematic uncertainties shown 
in Table 13 were estimated from variations to nominal fit results in numerous cross- 
checks. We consider variation of the efficiency shape, uncertainties in resonance 
parameters, presence of additional resonances etc. We note that with higher statis- 
tics this mode may offer one of the best measurements of the parameters of the 
ao(980). 



Table 3. Preliminary fit results for Xcl tt+tt rj Dalitz plot analysis. The 
uncertainties are statistical and systematic. 



Contribution 


Amplitude 


Phase (°) 


Fit Fraction {%) 


ao(980)±7rT 


1 





56.2 ±3.6 ±1.4 


/2(1270)r, 


0.186 ±0.017 ±0.003 


-118 ±10 ±4 


35.1 ± 2.9 ± 1.8 


arj 


0.68 ±0.07 ±0.05 


-85 ±18 ±15 


21.7 ±3.3 ±0.5 



The Dalitz plot for Xci K^K~Tr'^ decay and its projections are shown in 
Fig. 01 and for Xci ^ K^K^n^ in Fig. El We do a combined Dalitz plot analysis 
to these modes taking into account isospin symmetry on amplitudes and phases: 

a^.+ = aK*- — aii*a = a-j^ = a^,, 0^*+ = 4'K'- — <I>K'0 — 'kx^ = <pK*, 

aa(980)+ — aa(980)- — a.a{gsoY> = "a(980), ^ud 0q(98O)+ = 0a(98O)- = 0a(98O)O = 

</'o(980) ■ The common relative factor between two Dalitz plot amplitudes does not 
matter because of each Dalitz plot has an arbitrary normalization. The over- 
all amplitude normalization and one phase are arbitrary parameters and we set 
a-K' {892) = 1 and 0k* (892) = 0. 

The limited size of this sample, even in the combined Dalitz plot analysis, and 
the many possible contributing resonances leave us unable to draw clear conclu- 
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Fig. 4. Dalitz plot and projections on the two mass squared combinations for Xci ^ K^K tt'^ 
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Fig. 5. Dalitz plot and projections on the two mass squared combinations for Xci — > ''^^ ■ 

sions. Visual inspection shows clear contributions from iiT* (892)=*= iiT^, if *(892)°i4'^, 
i4:*(1430)±ifT^ K*{\^mfK%, aQ{^mfTP, and ao(980)±7r=F. It is not clear that the 
i^T* (1430) are or ifj, and many other K-n and KK resonances can possibly 
contribute. Our best fit preliminary result is presented in Table 01 showing statis- 



Table 4. Preliminary results of the combined fits to the Xci 
and Xci K^Ktt Dalitz plots. 



Contribution 


Amplitude 


Phase {°) 


Fit Fraction (%) 


ft:*(892)K 

K{{lim)K 
i^:* (1680) A' 
ao(980)7r 


1 

0.50 ±0.09 ±0.12 
5.3 ± 1.0 ±0.1 
2.3 ±0.5 ±0.5 
10.8 ± 1.2 ± 1.2 




-2 ±13 ±6 
77± 12 ± 16 
-38 ±12 ±12 
-112 ±12 ±3 


19.7 ±4.0 ±2.0 
18.0 ±6.6 ±3.2 
36.0 ±12.8 ±3.0 
11.2±5.4±2.7 
29.5 ±7.1 ±2.6 



tical and systematic errors. This fit has a good 32% probability of matching the 
data and agrees with fits done to the separate Dalitz plots not taking advantage of 
isospin symmetry. We also find that an alternative solution using the same set of 
contributions fits the data acceptably, but less well than the displayed result. 
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5. Summary 

The CLEO collaboration possesses data samples which can be used for analysis 
of D, Ds, Xcj, J/i' meson decays. Using roughly 10^ of D^D~ pairs produced 
in the e+e^ —t -0(3770) DD process we perform a Dalitz plot analysis of the 
D'^ — > 7r~7r"*'7r"'" decay. We confirm a dominant contribution from a low mass tt+tt" 
S-wave amplitude earlier observed by E791 and FOCUS experiments. Results of the 
_D+ — > 7r^7r+7r+ Dalitz plot analysis are summarized in Tabled 

We have searched for and studied selected three body hadronic decays of the Xco, 
Xcij Xc2 produced in radiative decays of the "0(25') in e^e~ collisions observed with 
the CLEO detector. The clean signal at least for one of Xcj's is found in eight final 
states: Tr+Tr^jy, tt+tt^?]', K^tt*^ , KgKir, rjK^K^, t^^pP, VPP^ AK^p. Our 
preliminary observations and branching fraction limits are summarized in Table |21 
In Xci — *■ 7r"'"7r~77 we have studied the resonant substructure using a Dalitz plot 
analysis, and our preliminary results are summarized in TableOl Similarly in Xci ^ 
KKtt we clearly see contributions from K*{892)K and ao(980)7r at roughly the 20% 
and 30% level respectively, as shown in Table 01 
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